Trapped surfaces and emergent curved space in the Bose-Hubbard model 
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A Bose-Hubbard model on a dynamical lattice was introduced in previous work as a spin system 
analogue of emergent geometry and gravity. Graphs with regions of high connectivity in the lattice 
were identified as candidate analogues of spacetime geometries that contain trapped surfaces. We 
carry out a detailed study of these systems and show explicitly that the highly connected subgraphs 
trap matter. We do this by solving the model in the limit of no back-reaction of the matter on 
the lattice, and for states with certain symmetries that are natural for our problem. We find that 
in this case the problem reduces to a one-dimensional Hubbard model on a lattice with variable 
vertex degree and multiple edges between the same two vertices. In addition, we obtain a (discrete) 
differential equation for the evolution of the probability density of particles which is closed in the 
classical regime. This is a wave equation in which the vertex degree is related to the local speed 
of propagation of probability. This allows an interpretation of the probability density of particles 
similar to that in analogue gravity systems: matter inside this analogue system sees a curved space- 
time. We verify our analytic results by numerical simulations. Finally, we analyze the dependence 
of localization on a gradual, rather than abrupt, fall-off of the vertex degree on the boundary of the 
highly connected region and find that matter is localized in and around that region. 
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I. INTRODUCTION 

Since the discovery of the Hawking and Unruh effects 
it has been clear that gravity is fundamentally different 
from the other forces. That a new thermodynamics has 
to be associated to black hole physics is a remarkable 
puzzle which physicists are slowly unveiling. In the last 
two decades, the possibility that gravity itself may have a 
thermodynamical origin has been explored from a variety 
of angles pUI] . A closely related idea is that gravity may 
be emergent, either the thermodynamics of a microscopic 
(quantum) theory [5-7 j, or emergent in the condensed 
matter sense [BUTT]. 

Motivated by the possibility that gravity may be emer- 
gent and the questions this raises, quantum graphity was 
introduced in as a method to study the emergence 

of geometry and gravity in the simplified but explicit and 
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tractable setting of a spin system. Graphity models are 
systems of quantum, dynamical graphs where the spin 
degrees of freedom correspond to dynamical adjacency: 
the edges of the graphs can be on (connected), off (dis- 
connected), or in a superposition of on and off. By inter- 
preting the graph as pre-geometry (the connectivity of 
the graph tells us who is neighboring whom), this pro- 
vides a way to treat pre-geometry as a spin system. A 
particular graphity model is given by such graph states 
evolving under a local Ising-type Hamiltonian. The goal 
is to investigate whether such a local dynamics can lead 
to the system exhibiting aspects of gravity in the ther- 
modynamical limit. In previous work we were able to 
derive features such as a variable speed of light from a 
non-geometric quantum system; we studied the quantum 
dynamics in a spin system as a precursor to a quantum 
theory of gravity ; we showed entanglement and thermal- 
ization between geometry and matter; and found that 
certain states appear to behave like geometries contain- 
ing trapped surfaces. 

These models can also be viewed as analogue mod- 
els for gravity. In the analogue gravity program vari- 
ous phenomena of general relativity (e.g., black holes or 
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cosmological geometries) are modeled by other physical 
systems, such as acoustics in a moving fluid, superfluid 
helium, or Bose-Einstein condensate; gravity waves in 
water; and propagation of electromagnetic waves in a di- 
electric medium [3J. The fluid analogues provide a very 
interesting concrete setting in which to study emergent 
gravitational properties, but leave open several impor- 
tant questions: 1) What is the role of the background 
fluid? 2) Can we go beyond analogues of kinematical 
only aspects of general relativity? 3) Lorentzian geom- 
etry: in analogue gravity systems, Lorentzian geometry 
"emerges" in the sense that it is the geometry seen by 
excitations inside the fluid (phonons). What happens 
if there is no background fluid, and can this work in a 
quantum system? 

Alternative analogue models, including ours, may at- 
tempt to answer some of these. For example, in our 
model there is no background fluid. There is, however, 
a background time in the Hamiltonian. Interestingly, as 
we will show in this paper, matter in our model sees an 
effective Lorentzian geometry, a phenomenon that can be 
traced to the locality of the Hamiltonian. 

In this article, we gain substantial insight into the 
model of |13j by studying two important questions: we 
explicitly confirm the argument of |13) that highly con- 
nected regions in the graph trap matter, and we relate 
the evolution of matter under a Hubbard Hamiltonian on 
these lattices to a wave equation on a curved geometry. 
We achieve this by restricting the model to a large class 
of states relevant for these questions. These are states 
in which the graph has certain symmetries (although the 
vertex connectivity and the number of edges linking two 
vertices can vary) and states of matter that respect these 
symmetries. 

We find that the Bose-Hubbard Hamiltonian with ho- 
mogeneous couplings on a lattice with varying vertex 
degree is equivalent to a Bose-Hubbard Hamiltonian on 
a regular-degree graph but with site-dependent effective 
couplings (similar to the behavior seen in [2]). This 
makes it possible to connect the coefficients in the Hamil- 
tonian to geometric properties of the graph. The picture 
which emerges from this analysis is the one of Fig. [T] 
the graph modifies the strength of the interaction in the 
Bose-Hubbard Hamiltonian, and this appears as a curved 
geometry to the propagating matter. 

It is important to note that the emergent curved space 
is a dynamical property of the system. The geometry that 
the particles propagating on the graph see depends on the 
dynamics of the particles and it is not just a property of 
the graph, and, in addition, the resulting motion of the 
particles will change the graph and so affect the geometry. 

The structure of the paper is as follows. In Section 
[n]we review the model introduced in [T3] and state the 
approximations we will use. This amounts to a Hub- 



bard model on a fixed but unusual lattice. In Section III 




J 



FIG. 1. The scheme representing the relations between the 
graph, the hopping energies fi t j of the Bose-Hubbard model 
and the emergent metric: the graph modifies the strength of 
the interaction in the Bose-Hubbard Hamiltonian, which in 
turn translates into a curved geometry (for the appropriate 
states) . 



WJ we study the states that were previously identified 
as possible trapped surfaces and confirm this property. 
These states have the symmetries defined in the previ- 
ous Section. In the single particle sector we can diag- 
onalize the Hamiltonian on the subspace of rotationally 
invariant states and show that the ground state belongs 
to that subspace. We then show that this ground state 
is protected by a gap which increases linearly with the 
size of the completely connected subgraph, so when that 
size is large, this area acts as a trap for matter. In the 
many-particle sector, we find that the ground state is a 
Bose-Einstein condensate of delocalized particles at the 
completely connected region, again with a large gap. We 
confirm these results numerically. In Section |VJ we re- 
late the one-dimensional Bose-Hubbard model with site- 
dependent coefficients to the wave equation in curved 
space for the particle probability on the lattice. In sec- 
we study numerically matter localization as a 
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we define foliated and rotationally invariant graphs, the 
symmetric graph states we will consider, and delocalized 
states, our similarly symmetric matter states. In Section 



tion 

function of a fall-off parameter of vertex degree on the 
boundary of the highly connected region. We find that 
the particles are localized inside the trapped surface as 
before, but also on a small region around it. Conclusions 
follow. 



II. THE MODEL 

In this section we will introduce the model, starting 
with the basic idea of describing a dynamical lattice using 
quantum degrees of freedom on the edges of a complete 
graph. 

We start with a quantum mechanical description of 
a universe of N v elementary systems, given by the set 
{Jii, Hi} of the Hilbert spaces JC, and Hamiltonians Hi of 
the systems i = 1, N. This presumes it makes sense to 
talk of the time evolution of an observable with support 
on Hi without any reference to a spatial coordinate for 
i. 

We choose J£j to be the Hilbert space of a harmonic 
oscillator. We denote its creation and destruction op- 
erators by b\ and bi respectively, satisfying the usual 
bosonic commutators. Our N v physical systems then are 
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N v bosonic modes and the total Hilbert space of such 
modes is given by 

N v 

~N-boso?is (1) 
i=l 

If the harmonic oscillators are not interacting, the total 
Hamiltonian is trivial: 

N v JV„ 

H v =J^H i = -J^ f xb\b i . (2) 

i=i i=i 

If, instead, the harmonic oscillators are interacting, we 
need to specify who is interacting with whom. Let us 
call I the set of the pairs of oscillators e = (i, j) that are 
interacting. Then the Hamiltonian reads as 

i eel 

where h e is a Hermitian operator on Hi®Hj representing 
the interaction between the system i and the system j. 

Spacetime geometry determines the set of relations 
among physical systems, namely, interactions are only 
possible between neighbors. In a discrete setup, we can 
introduce a primitive notion of spacetime geometry via 
the adjacency matrix A, the N v x N v symmetric matrix 
defined as 

^ _ f 1 if i and j are adjacent 
lj ~ \ otherwise. ^ ' 

A defines a graph on N v nodes, with an edge between 
nodes i and j for every 1 entry in the matrix. Now note 
that all possible graphs on N v vertices are subgraphs of 
Kn, the complete graph on N v vertices (with adjacency 
matrix Ay = 1 for every i,j). In our model, the N V (N V — 
l)/2 edges of Kn correspond to (unordered) pairs e = 
of harmonic oscillators. To every such pair e (an 
edge of Km) we associate a Hilbert space J{ e — C 2 of a 
spin 1/2. The total Hilbert space for the graph edges is 
then 

N v {N v -l)/2 

^graph = (^) M e . (5) 
e=l 

The basis in ^K grap h is chosen so that to every subgraph 
G of Kn corresponds a basis element in ^K gra ph'- the 
basis element \e\ . . . ejv„(jv„-i)/2) = \G) corresponds to 
the graph G that has an edge e s present for every e s = 1. 
For every edge (i, j), the corresponding SU(2) generators 
will be denoted as S l — l/2er l where a z are the Pauli 
matrices. 

The total Hilbert space of the theory is 

^ ^-bosons ® ^ graph: (6) 

and therefore a basis state in H has the form 

= \rii,...,n N „) (g) \ei,...,e N v t.N v -i) ). (8) 



The first factor tells us how many bosons there are at 
every site i. The second factor tells us which pairs 
interact, that is, the structure of interactions is now pro- 
moted to a quantum degree of freedom. Unlike the model 
presented in |13j where hard core bosons were considered, 
the interaction among the bosons in the present paper is 
an on-site potential described by ub\bi(b\bi — l)/2 where 
u is the interaction energy between two particles. 

In general, a quantum state in the full model describes 
a system in a generic superposition of energies of the har- 
monic oscillators, and of interaction terms among them. 
A state can be a quantum superposition of "interactions" . 
For instance, consider the systems i and j in the state 

I*,) = »M)i_k, (9) 
v2 

This state describes the system in which there is a parti- 
cle in i and no particle in j, and a superposition of link 
and no link between i and j. 

A simple but interesting matter interaction term is the 
one that describes the physical process in which a quan- 
tum in the oscillator i is destroyed and one in the os- 
cillator j is created. This dynamical process is possible 
when there is an edge between i and j. Such dynamics 
is described by a Hamiltonian of the form 

H hop = -E hop % ® (blbj + b^), (10) 

where 

p v = = i^icu) = Q - ( n ) 

is the projector on the state such that the edge (i, j) is 
present and the spin operators are defined as ^uj) ~ 
|l)(0|(ij) and Suj) = |0) ( 1 1 - With this Hamiltonian, 
the state defined in Eq. ([£]) can be interpreted as the quan- 
tum superposition of a particle that may hop or not from 
one site to another. 

It is possible to design such systems in the laboratory. 
For instance, one can use arrays of Josephson junctions 
whose interaction is mediated by a quantum dot with two 
levels. 

We note that it is the dynamics of the particles de- 
scribed by Hhop that gives to the degree of freedom |e) 
the meaning of geometry. The geometry at a given in- 
stance is given by the set of relations describing the dy- 
namical potentiality for hopping. Two vertices j, k can 
be "empty" , that is, the oscillators j, k are in the ground 
state, but they can still serve to allow a particle to hop 
from i to j, then to k, then to I. We read out the structure 
of the graph from the interactions, not from the mutual 
positions of particles. 

In addition, -Hhop tells us that it takes a finite amount 
of time to go from i to j. For instance, in a one- 
dimensional graph (a chain) it takes a finite amount of 
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FIG. 2. Diagramatic representation of the exchange and hop- 
ping terms of the model. The exchange term destroys a link 
and creates two particles on the vertices of that link; the par- 
ticles hop; two particles can combine to form a new link. 

time for a particle to go from one end of the chain to 
another (modulo exponential decaying terms) . As shown 
in [TS] , this results in a "spacetime" picture with a finite 
light-cone structure. 

The hopping amplitude is given by t, and therefore all 
the bosons have the same speed. Note that, for a larger 
Hilbert space on the links, we can have different speeds 
for the bosons. 

Of course, we need a Hamiltonian also for the spatial 
degrees of freedom alone. The simplest choice is to assign 
an energy to every edge: 

Hunk = -uJ2^i,jy ( 12 ) 

Finally, a central motivation for the model [13] was to 
have space and matter interact unitarily. The term 

H ex = k £ ® (b\b]) Q + S+j) ® W 3 ) ' ( 13 ) 

(id) 

introduced in |13| . destroys an edge and creates Q 
quanta at i and Q quanta at j, or, vice-versa, destroys Q 
quanta at i and Q quanta at j to convert them into an 
edge (see Fig. [2| . 

As mentioned above, the long-term ambition of these 
models is to find a quantum Hamiltonian that is a 
spin system analogue of gravity. In this spirit, matter- 
geometry interaction is desirable as it is a central feature 
of general relativity. The above dynamics can be consid- 
ered as a very simple first step in that direction. 

This was the final step that brings us to the total 
Hamiltonian of the model which is 

H = Hunk + H v + H ex + Hbop- (14) 

This is the original model proposed in [T3] . It is a Bose- 
Hubbard model on a dynamical lattice, where the pro- 
jector on the modified hopping term means that the dy- 
namics gives the lattice the meaning of space: space is 
where matter is allowed to go. 

In the present work, we study the model for a par- 
ticular class of graphs that have been conjectured to be 
analogues of trapped surfaces. We are interested in the 
approximation fc <C t, which can be seen as the equiv- 
alent of ignoring the backreaction of the matter on the 
geometry. For simplicity, we will set U = k = 0, meaning 
that 

H = H V + H hop . (15) 



In this case, the total number of particles on the graph is 
a conserved charge. H v and Hu n k s are constants on fixed 
graphs with fixed number of particles. The Hamiltonian 
is the ordinary Bose-Hubbard model on a fixed graph, 
but that graph can be unusual, with sites of varying con- 
nectivity and with more than one edge connecting two 
sites. Our aim will then be to study the effective geome- 
try that matter in this model sees. 

Even on a fixed lattice, the Hubbard model is diffi- 
cult to analyze, with few results in higher dimensions. It 
would seem that our problem, propagation on a lattice 
with connectivity which varies from site to site is also 
very difficult. Fortunately, it turns out that for our pur- 
poses it is sufficient to restrict attention to lattices with 
certain symmetries. 

III. FOLIATED AND ROTATION ALLY 
INVARIANT GRAPHS 

In this Section we define graphs with two particular 
properties which we call foliation and rotational invari- 
ance. These properties will allow us to greatly simplify 
the calculations that follow without loss of generality. We 
will see that the problem of finding the ground state of 
hopping Hamiltonians, on graphs with these properties 
can be simplified to the solution of the one-dimensional 
Bose-Hubbard model. 

A. Foliated graphs 

A foliated graph is a graph that can be decomposed 
into a set of subgraphs connected by edges in a successive 
way. More precisely, let gi be a labeling of subgraphs of a 
graph G and Ei a labeling of the set of edges connecting 
the sets gi, such that U{(<fc U Ei) = G. Then, 

Definition 1 A graph G is foliated if it can be decom- 
posed in several disjoint subgraphs gi with the following 
properties: 

1. All the subgraphs gi are degree regular. 

2. All the edges of a subgraph gi connect a vertex in 
gi to a vertex in gi-i or 

3. The number of edges connecting a vertex in gi 
(gi+i) to vertices in gi+i (gi) is the same for every 
vertex of gi (gi+i)- This number is called relative 
degree and is represented by dj -j+i (di+xj). 

Notice that the name, foliated, comes from the fact that 
these graphs can be decomposed into subgraphs such as 
any foliated structure can be separated into thin layers. 
Examples of foliated graphs are presented in Figs. [3] and 

m 

The number of edges that connect two consecutive sub- 
graphs gi and g i+1 is given by 

d l . i+1 N i = d i+1 . t N. l+ i , (16) 
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FIG. 3. A 2-d lattice that can be foliated. 




FIG. 4. A non-planar lattice that can be 1-foliated. 



where iV, is the number of vertices of gi. 

In order to discuss the properties of the hopping Hamil- 
tonians defined by foliated graphs, let us introduce a nat- 
ural coordinates (i, m) to specify a vertex of a foliated 
graph. The coordinate i specifies the subgraph gi and m 
refers to the particular vertex in gi. The hopping Hamil- 
tonian of a foliated graph in terms of these coordinates 
can be written as 

R-l Ni-l 

^° 1 = -E E A m m 4,mKm>+h.C. 
2—0 m,m ; — 
R-l Ni-l N i + 1 -1 
"EE E B{ ±4,mh + l,m'+h.C, (17) 
j=0 m=0 m'=0 



where b\ m (bi im ) is the creation (annihilation) operator of 
a particle at the vertex (i, to), R is the number of layers 
9ii ^-mrn' ^ s the adjacency matrix of the subgraph gi, and 

(i) 

"mm 1 stands for the edges between the vertices of gi and 
9i+i- 

In this model, the delocalized states on the subgraphs 
gi are particularly interesting. These are states of a par- 
ticle which is completely and uniformly spread over the 
graph g i: defined as follows: 



Definition 2 The delocalized state, in gi is defined 
by 



Ni — l 



^E 

/7V7 ^ 



(18) 



The delocalized state |i) is an eigenstate of the hopping 
Hamiltonian defined from the degree-regular subgraph gi, 



H, 



« 
fol 



JVj-1 

E 

m.m'—O 



A 



ram' i.in t ; rn ' 



+ h.c. 



(19) 



More explicitly, 




-diiE, 



hop | 



(20) 

where we have used the degree-regularity of the graph, 

= E m ^™' ^ 0r au n - ^ 0r example, when the graph 
gi is a chain, the completely delocalized state is the 
ground state of the system and has energy —2, with 2 
is the degree of the vertices of a 1-dimensional chain. 

The subspace spanned by the delocalized states \i) is 
an eigenspace of the system, since the projector onto this 
subspace, 



Pa 



R-l 

Ei 

i=0 



(21) 



commutes with the Hamiltonian (17 1, [Hf a \, P<j] = 



Therefore, the time evolution of any superposition of de- 
localized states lies always in this subspace (it is a super- 
position of delocalized states). This allows us to define 
an effective Hamiltonian for the delocalized states, 

H cS := P d H io iP d (22) 

R-l 

= E (A** + + 1^ + !D + wl<X<l) . 

i=0 

where /i, = (i\H iol \i) = -E hop d t j and 



fi,i 



1+1 



(i|flfoi|i + 1) = -^ho P di,»+i 




(23) 



Note that the Hamiltonian ( 22 1 is a one-dimensional 



m=0 



Bose-Hubbard Hamiltonian with chemical potential \ii 
and tunneling coefficient /j^+i. 

The mass term (or chemical potential) /i, is fattened 
by the edges connecting the nodes within the subgraph gi 
of the foliated graph. This is one the main results of 
the paper. This behavior resembles a scalar field not- 
minimally coupled to (classical) gravity, where the mass 
of the particle is multiplied by a curvature factor. In our 
case, the role of the curvature is played by —Eh op di t i (see 
SecjV}. 

The extension to higher dimensions is straightfor- 
ward. It requires the extension of one-dimensional fo- 
liated graphs to graphs which can be foliated in multiple 
directions, thus resembling an ordinary lattice, but with 
multiple links between pairs of sites. The coefficients 
fkk-i w iU depend on the direction of the foliation that 
the particle is hopping to. 



We have found an eigenspace of the foliated Hamilto- 
nians for which the effective Hamiltonian is particularly 
simple. However, this eigenspace will not, in general, con- 
tain relevant eigenstates of the Hamiltonian such as the 
ground state. In order to ensure this, we require another 
symmetry. 



B. Rotationally invariant graphs 

Let us present next our definition of the rotationally 
invariant graphs. 

Definition 3 A graph G is called N -rotationally invari- 
ant if there exists an embedding of G to the plane that is 
invariant by rotations of an angle 2ir/N . 

With the exception of those graphs that have a vertex 
in the center, N- rotationally invariant graphs allow for a 
labeling of their vertices (n,9) € N x Zjv such that their 
adjacency matrix A^g^ ^ n igi) only depends on n, n 1 and 
8 — 6'. We can make use of these coordinates (n, 9) in 
order to write the Hamiltonian defined by a rotationally 
invariant graph 



N-l 



H Iot = - ^2 ^2 A nn ,b\ e b nlf) + h.< 



0=0 n,ri 
N-l N-l 



-EEE*^w + k - ( 24 ) 

9=0 ip=l n,n' 

where b* n g (b n .e) is the creation (annihilation) operator 
at the vertex (n,6), A nn i is the adjacency matrix of the 

graph of any angular sector and i?i „/ is the adjacency 
matrix of two angular sectors at an angular distance ip 
in units of 2tt/N. 

Let us introduce the rotation operator L defined by 



Lb n fi = b ry + iL 



Lb 



(25) 



The effect of the operator L is particularly easy to un- 
derstand in the single particle case: 



Lin, 



Lbl e \0} = bl d+1 L\0) = \n,e + l) 



(26) 



where we have assumed that the vacuum is invariant un- 
der a rotation L|0) = |0). 

Note that L is unitary and its application N times gives 
the identity, L N = 1. This implies that its eigenvalues 
are integer multiples of 2tt/N. 

Another interesting property of L is that commutes 
both with the rotationally invariant Hamiltonians and 
with the number operator N p , 




FIG. 5. A foliated graph which is not rotationally invariant. 



Therefore H mtl N p , and L form a complete set of com- 
muting observables and the Hamiltonian is diagonal in 
blocks of constant L and N p . This allows us to simplify 
the problem of finding the ground state and the first ex- 
cited states of the system, a very useful fact we will use 
in the next section. 

Note that rotational invariance and foliability are not 
equivalent. There are graphs which can be foliated and 
are not rotationally invariant and vice- versa. An example 
is given in Fig. [5] 

We have seen that both for the cases of foliated 
graphs and the rotationally invariant graphs, the sub- 
space spanned by the completely delocalized states is an 
eigenspace of the system. In the foliated model, ef- 
fective dynamics in this eigenspace are given by a one- 
dimensional Bose-Hubbard Hamiltonian. In the rota- 
tional invariant case, this eigenspace contains the ground 
state of the system. Thus, the ground state of graphs 
with both symmetries lies in the subspace of completely 
delocalized states and the computational effort to find it 
is equivalent to the solution of a one-dimensional Bose- 
Hubbard model. This allows us to analyze a complicated 
model using the approximation of a one-dimensional spin 
chain, with obvious advantages, especially for numerical 
work. 



IV. 



REGIONS OF HIGH CONNECTIVITY AS 
TRAPPED SURFACES 



In [T3], we observed that the Hamiltonian (15) caused 



[H mt , L] = [N p , L] = [H mU N p ] = . 



(27) 



trapping of matter in regions of higher connectivity. The 
basic mechanism is the following: consider a graph con- 
sisting of two set of nodes, A and B, separated by a set 
of points C on the boundary. Let the vertices in A be 
of much higher degree than the vertices in B, dA 3> ds 
(see Fig. [6| . If the number of edges departing from the 
set C and going to the set A is much higher than the 
number of edges going from C to B, then the hopping 
particles have a high probability of being "trapped" into 
the region A. 

Our task in this paper is to make this heuristic ar- 
gument precise and determine whether these high con- 
nectivity configurations are spin-system analogues of 
trapped surfaces. We will do this by studying specific 
states that are graphs with symmetries that contain a 
core (trap) of N nodes. Fig. [7] is an example of such a 
graph. Such states are I-foliated graphs and we will be 



FIG. 6. A region of higher connectivity in a regular graph. 



FIG. 7. The X N graph. 



able to use the properties we discovered above. 



A. Classical trapping 

In order to gain some intuition on the trapping, let us 
consider the classical analogue of the problem. In the 
classical setup, a particle has a well-defined position in 
some site of A. At each time step of length h/Eh p, 
the particle randomly jumps to another site of the graph 
connected to its current site by an edge. This process is 
successively repeated until the particle escapes from A. 
How much time is required (on average) for the particle 
to escape from the highly connected region? 

In this simple model, the probability that the particle 
jumps to a site outside A is 



Pc 



N-l + N c: 



(28) 



where N oxt is the number of links that connect a site in 
A with the environment, and N is the total number of 
sites in A. The evaporation time is then given by 



to 



1 



N 



-^hop Pesc 



-^hop -^ext 



(29) 



where we have assumed that N ^> N cxt - In the large N 
limit, the particle is trapped in A. 



B. Quantum case: the %n configuration 

In this section we study the spectrum and the con- 
figuration of the ground state and the finite temperature 
states of the hopping Hamiltonian on the graph of Figure 
[7] which we will call the 3Cjv graph. This is a 1-foliated 
graph with a completely connected core. We will show 
that the model on %m can be solved analytically in the 
thermodynamic limit. 

The position of a vertex of the %n graph can be speci- 
fied by means of the integer spherical coordinates r and 9, 



with ranges < r < R- 1 and < 9 < N- 1. Then, the 
quantum state of a particle with a well-defined position 
in the graph can be written as 



\r,e) = bljo) 



(30) 



where |0) is the vacuum state and b r 6 the corresponding 
creation operator. Using these coordinates, the hopping 
Hamiltonian defined by the 3Cj\r graph becomes 

N-l E-i 

H a =J2J2 ( 6 H-iA,A + b i,e+A,e + h - c ) ( 31 ) 

9=0 r=0 

+ b ke' b o,e, 

\8-8'\>2 

where h.c. is the Hermitian conjugate. Note that the 
second term in Eq. (31) corresponds to hopping in the 



completely connected region, while the first sum is hop- 
ping outside that core. Our question is how the intro- 
duction of this completely connected region (the second 
sum) changes the spectrum and the eigenstates of the 
system. 

Single particle case. Let us first work out the single 
particle sector of the Hamiltonian. 

In order to determine the eigenstates and eigenvalues 
of Hq, we will write the Hamiltonian in the eigenbasis 
of the rotation operator L defined in equation (25). The 
eigenstates of L read 



N-l 



with eigenvalues 



?=0 



L\r,i) =e 1 ^ i \r,£), 



(32) 



(33) 



where t — 0,1, . . . , N — 1. The Hamiltonian is diagonal 
in blocks of constant I and can be written as 



N-l 



(34) 



where Hi — PiHPi are the projections onto the 
eigenspaces of L with the projectors Pi defined as 



R-l 



(35) 



r=0 



fnserting Eq. (35) into the definition of He, we get 

'2tt 



H e = -N5 m \Q)(Q\ - 2 cos 



N 



(|r + l)<r| + |r)<r + l|), 



(36) 



where 0(-) is the Heaviside step function, introduced to 
make the second term of the right hand side vanish when 
r = 0. 

Note that Hi is translationally invariant in the limit 
R . — >• oo for £ > 0. Therefore, it can be analytically 
diagonalized by the discrete Fourier transform 



\k r ,£) 

with k r = 0, . . 
H t \k r ,i)=2 



1 R-l 

1 V T : 2 



(37) 



R 



1 



r=0 

1 and 1 < £ < N 
'2vr 



cos 



R 



cos 



1. We find 



(38) 

up to 1/R corrections that vanish in the infinite limit. 

Let us now consider the subspace £ = of the rota- 
tionally invariant states. Because the Hamiltonian com- 
mutes with L, the ground state \GS) of the system must 
be invariant under its action: L\GS) = \GS). Therefore, 
\GS) belongs to this subspace. On it, we can explicitly 
construct the matrix for the Hamiltonian He = o'- 



I N 1 



Ho 



-Eh 



hop 



l 

o 

V o 







f 

1 2/ 



(39) 



It is a tridiagonal matrix with characteristic polynomial 
PrW — dct(H — XI) which can be written in a recursive 
way as 

Po(A) = 1, 

Pi(X) = -N, (40) 

Pn(A) = -(2 - X)p n -i(X) - p„_ 2 (A) . 

Note that because of the recursive relation, it is not clear 
whether the other eigenvalues apart from the first one 
depend on N or not. 

It is easy to see that, if we rescale Shop = Eh op /N and 
take N — ¥ oo, the only element left in the matrix is the 
element associated to the |0) state. Thus, at N — > oo, the 



ground state becomes |r = 0, £ = 0) and the gap between 
it and the first excited state scales as N. In the ther- 
modynamic limit, the ground state of the single particle 
sector corresponds to a particle completely dclocalized in 
the complete graph, 



N-l 



N 



lim \GS) = |r = 0,* = 0) 



9=0 



0) ■ (41) 



The rest of eigenvectors of the subspace £ — are or- 
thogonal to \GS) = \r — 0, £ = 0), and therefore lie in the 
subspace spanned by \r,£ = 0) with r > 1. The Hamil- 
tonian can be analytically diagonalized in this subspace 
by the same transformation used in Eq. ( 37 ) : 



R-l 



(42) 



withfc = 0, 2. 

In conclusion, we have seen that the %n model has a 
unique ground state which is protected by a gap which 
increases linearly with the size N of the completely con- 
nected region. The rest of the eigenvalues form an energy 
band which is the same as if we had the %n graph with- 
out the completely connected region, 



Ekt = 2(1- cos 



2tt. 
R 



k r I — cos 



2tt 
— i 
N 



(43) 



These results are numerically confirmed by Figs. [8] and 
[9] In Fig. [8j the energies of the ground state, the first 
excited state and the state with maximum energy are 
plotted against the size of the completely connected re- 
gion. In Fig. [9j we have plotted the fidelity between the 
ground state of the system and the completely delocal- 
ized state in the completely connected region vs N. 

Multi-particle case. Let us next analyze what happens 
when there are several particles in the %n configuration, 
interacting with an on-site potential. 

The Hamiltonian of the multi-particle 5Cjv model can 
be decomposed in its one-body and two-body parts, 



H = H Q + V 



(44) 



where Hq is defined in Eq. (31) and V is an on-site in 
teraction among the particles: 

R-l N-l 



V = u 



£ J2 h le h le h reKe , 

r=0 6=0 



(45) 



with u the energy penalty for two particles in the same 
site. 

Because the interaction V commutes withj;he_number 
operator iV and the rotation transformation L, H, L and 
TV form a complete set of commuting observables. It is 
convenient to write Hq and V in terms of the creation 



9 



20 



a 



-20 



-40 



-60 




10 



20 



:!0 



40 



50 



FIG. 8. Plots of the energy Eo of the ground state, the energy 
Ei of the first excited state, and -E ma x, the energy of the 
maximum energy state, against the size N of the completely 
connected region of the %n graph, in the single particle sector. 
This plot has been realized for the full model with R — 30. 
Note that the gap E± — Eo increases linearly with N. 
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FIG. 9. Fidelity between the completely delocalized state 
in the completely connected region |r = 0, £ = 0) and the 
ground state against the size N of that region. The completely 
delocalized state |r = 0,1 = 0} becomes the ground state of 
the system for large N . 



the Hamiltonian is 



H = -N V 1 V; 



gs'/gs 



5>< 



(49) 



and the interaction reads 



kl 



+ u vLvtsVgsrigs , 



(50) 



where the sum runs over all possible k\ and ti for i = 
1, . . . , 4 and the Kronecker delta shows that the interac- 
tion conserves the quantum numbers I and k. Note that 
the single-particle ground state is completely decoupled 
from the states of the energy band. Therefore, the state 
of n particles, 

\GS) = ( V l) n \0), (51) 

is an eigenstate of the system with energy ~nN + un(n — 
l)/2. Furthermore, as the interaction is a positive oper- 
ator, we can ensure that \GS) is the ground state of the 
system as long as n < N. 

Thus, the ground state of the many-body problem is a 
Bose-Einstein condensate of delocalized particles at the 
completely connected region. The large gap and the 
features of the on-site interaction make this condensate 
robust at finite temperatures and against adding inter- 
acting particles. Thus, the completely connected region 
of the %n model can be considered a trapped surface. 

However, it is not clear whether every completely con- 
nected region is a trapped surface. In particular, we 
would like to see what happens when the connectivity 
does not change as abruptly as in the %n system. For 



this reason, in Section VI we will parametrize the fall-off 
of the parameters and study the localization of particles. 
Moreover, we would like to relate this fall-off to an effec- 
tive curved space-time geometry, as in an analogue grav- 
ity system. In order to do this, then, a correspondence 
between the connectivity of the graph and the curvature 
of the space-time is required. Establishing this relation 
is our task for the next section. 



V. CORRESPONDENCE BETWEEN GRAPH 
CONNECTIVITY AND CURVED GEOMETRY 



and annihilation operators 

JV-l R-l 



VkO 



Vgs 



= EE ei * ferei ^M< w ^ ( 46 ) 

(47) 
(48) 




As it has been previously shown, the one-body term in 



In this section we establish a relation between the con- 
nectivity of a graph and the curvature of its continu- 
ous analogue geometry. In order to do this, we restrict 
the time-dependent Schrodinger equation to the mani- 
fold formed by the classical states, that is, single-particle 
states with a well-defined but unknown position. The 
equation of motion obtained corresponds to the equation 
of propagation of light in inhomogeneous media, simi- 
larly to black hole analogue systems. Once we have such 
a wave equation, we can extract the corresponding met- 
ric. In the second part of the section, the dispersion 
relation and the continuous limit are discussed for the 
transitionally invariant case. 
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A. Restriction of the time-dependent Schrodinger 
equation to the set of classical states 

Since we want to study the dynamics of a single particle 
on a fixed graph, it is only necessary to consider the sin- 
gle particle sector. The one dimensional Bose-Hubbard 
model for a single particle reads, 

Af-l 

Ho=Y, fn,n+i(\n)(n + l\ + \n + l)(n\)+J2^n\n)(n\, 

n—0 n 

(52) 

where f n ,n+i are the tunneling coefficients between sites 
n and n + 1, p, n is the chemical potential at the site n, 
and M is the size of the lattice. 

In this setup, let us introduce the convex set of classical 
states M.C-, parameterized as 



M-l 



(53) 



where 9 n is the probability of finding the particle at the 
site n. The states in Aic are classical because the uncer- 
tainty in the position is classical, that is, they represent 
a particle with an unknown but well-defined position. 

The aim of this section is to restrict the Schrodinger 
equation of the whole convex set of density matrices to 
the convex set Aic and obtain the effective equations of 
motion for the classical states. In order to do this we 
will follow the same procedure as in Ref. [H]. We will 
approximate the time evolution generated by Hq without 
ever leaving the convex set Mc- This procedure consists 
basically of two steps: 

1. Time evolution. Insert the initial state p(t) € Mc 
into the time-dependent Schrodinger equation to 
get its evolution p[t + At) after a short time At. 

2. Restriction to Mc- Find the state in Mc that best 
approximates the evolved state p(t + At). 

If we take the infinitesimal limit At — > of the previous 
steps we are going to get a differential equation for the 
field * n (i). 

Let us mention that the time-dependent Gross- 
Pitaevskii equation can also be derived by this method, 
that is, by restricting the Schrodinger equation to the 
manifold defined by states parametrized by \<p) = 
exp (J dx(p(x)b x ) |0), where b x is the field operator that 
creates a particle at position x. In Ref. [TB]> the second 
step is something that must be done in order to restrict 
the equations of motion to the desired manifold. Never- 
theless, in our case, decoherence can give a physical in- 
terpretation to this step. Since our particle is under the 
effect of a noisy environment, its density matrix is going 
to be constantly dephased by the interaction between the 
particle and its reservoir. 

If we insert p(t) in the time-dependent Schrodinger 
equation we obtain: 

d f p(t) = J2dt*n(t)\n)(n\ = i[H ,p(t)] , (54) 



where, in general, the right hand side cannot be written 
in terms of the left hand side, and therefore the previous 
equation cannot be fulfilled. Note that the Hamiltonian 
pushes p(t + At) out of the Mc- 

Thus, the best approximation to p(t+At) in the convex 
set Mc is obtained by minimizing the distance 



J2dt^n(t)\n)(n\-i\H ,p(t) 



(55) 



For the Hamiltonian (52), the commutator [Hq, pit)] 
reads 



[Ho At)] 



(56) 



= - *«) (\ n )( n + M + \n+ l)(n\) , 

n 

and therefore the state in Mc that best approximates 
p(t + At) fulfills 



d t V n (t)=0. 



(57) 



This forces us to consider the second order in At, and 
therefore the effective equation of motion for ^> n (t) is 
going to be a second order differential equation in time. 
We then have 



d? P (t) = j2 d ^^\ 



n)(n\ = — 



1 



H ,[H ,p(t)] 



The dephased state in the position eigenbasis that best 
approximates p(t + At) can be easily determined by com- 
puting the double commutator of the previous equation. 
It obeys the evolution 



(58) 



+ /^-l,„*n-l W + /n,n+l*n-lW. 

In order to rewrite the previous expression in a nicer way, 
let us add and subtract the quantity fn-i, n i^ n+i{t) + 
^ n (t)) in the right-hand side, obtaining 

^d^ n (t) =/*_ 1(B (* n+ i(t) + tf n _!(t) - 2* n (f)) 

+ (/"+!,„- fn-1, n) - *»(*)) ■ 

This equation becomes a wave equation in the continuum, 



d^(x,t) - d x (c 2 (x)d x ^>(x,t)) = , 



where 



1 

c(x) 



h 2 



(59) 



(60) 



and ty(x,t) and f(x) are the continuous limit functions 
of ^ n (t) and fn,n-i respectively. 
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We have shown that the equation of motion of the 
Bose-Hubbard model, restricted to the convex set formed 
by classical states, is the wave equation. This is a signif- 
icant result that establishes a relation between the cou- 
pling constants of the Bose-Hubbard model and the speed 
of propagation of the fields t). Equation (59) has 



the same form as the equation for propagation of light in 
media with a space-dependent refraction index, as is also 
the case in black hole analogue systems [17] . 

A few comments are in order. In equation ( 59 1 , the 
constant of the speed of light is quantized. It is propor- 
tional to the inverse of the number of links between the 
nodes in this simplified model. This constant depends on 
the hopping coupling constant of the hamiltonian, Shop- 
Finally, even though equation ( 59 ) is for a scalar quantity, 



the analogy with a Klein-Gordon field can not be pushed 
too far. The equation refers to the 1-particle density 
from which we started from, ^ n (t) = (n\p(t)\n) , and so 



is completely classical. A generalization of equation ( 59 ) 



to the case of many distinguishable interacting particles 
is given in detail in the Appendix, where we analyze also 
the effect of a local and non-local interaction. 



B. Dispersion relation and continuum limit 

Let us consider in more detail the translational invari- 
ant case in which / n _i, n = / and p n = fi for all n. In 
this case, the continuous wave equation ( 59 ) becomes 



d 2 V(x, t) - c 2 d^(x, t) = , 



(61) 



where c is the speed of propagation. 

In order to understand the continuum limit, we evalu- 
ate the dispersion relation for the propagation of proba- 
bility in the translationally invariant case. Let us intro- 
duce a discrete Fourier transform in the spatial coordi- 
nate and a continuous Fourier transform in the temporal 
coordinate, given by 



*n(<) - 



1 



M-l 

E 

fe=0 



*fc(t)e" 



(62) 



and *fc(t) = Ae iUkt + Be~ iwkt . After a straightforward 
calculation, we find that the relation between u>h and k 
is given by 



WfcC : 



V2 



'27T 

cos [ — k 
M 



(63) 



Note that this dispersion relation is different from that of 
the quantum excitations. In order to see this, we Fourier 
transform the Hamiltonian in the translational invariant 
modes. We define the field momentum vE'fc(t) as ^k(t) — 
(k\p(t)\k), where |jfc) = J2 m e i27rkm/M /VJ\f\m) are the 
eigenstates of the translationally invariant Hamiltonian 
H \k) = huj"\k), with 



huff = 2(1 



cos 



/27rfc\ 



(64) 



Therefore, 

*fc(f) = Y,j 2 " Hn ~ m) ' M P",™(t), (65) 
n,m 

and, for the classical states p nm {€) = 8 n ,m^nif), 

1 M_1 1 
**(*) = M E *»(*) = m- (68) 



Note now that Eq. (64 1 differs from (63). In fact, at the 
quantum level, the continuum limit gives the ordinary 
galilean invariance, while at the quantum level, the con- 
tinuum limit gives excitations (time-evolving probability 
densities) which are Lorentz invariant. 

However, the continuum limit is tricky because we have 
continuous time and discrete space, and no spatial scale 
to send to zero with the inverse of the number of lat- 
tice points (the graph is not embedded in any geome- 
try) . Thus the continuum limit is in the time scale of the 
modes. If we rescale ujk — > &k/M (or equivalently c), we 
find that 




2, 
M 



(67) 



and, therefore, 



lim Uk ~ 2tt 



k 



M- 



That is, only modes that are slow with respect to the time 
scale set by c see the continuum. Note that by rescaling 
the speed of propagation c, the continuum limit can be 
obtained by a double scaling limit, Shop — > Eh op /M and 
M — > oo for lattice size M. In this limit, the probability 
density has a Lorentz invariant dispersion relation. An- 
other interesting rescaling which gives Lorentz invariant 
dispersion relations is the rescaling of h. 

To conclude the discussion, let us note that the speed 
of propagation of the probability sets the timescale of the 
interaction: if there are P points with constant speed of 
light given by c, then the timescale of the propagation 
of the interaction in the classical limit in that region is 
given by t dyn = P/ c. 

An embedding of the graph in space, on the other hand, 
would mean requiring that ^ n (t) depends on a point in 
space and a length scale a (the lattice spacing) : ^ n (t) = 
^(60 + na,t), and that the coupling constant scales as 
Shop — > Sho P a. In this case the continuum limit a — > 
gives the ordinary dispersion relation. 



VI. MODEL WITH PARAMETRIZED FALL-OFF 
OF CONNECTIVITY 

In this section, we study a graph with a trapped sur- 
face (completely connected region) whose boundary is 
extended and the connectivity of its nodes decreases 
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gradually towards the outer edge of the boundary. We 
parametrize the fall-off of the parameters of the model 
and study the localization of a one-particle state. We 
use the techniques developed in Sees. |III| and |TV| to find 
the ground state of the system and check its robustness p, - 



against finite temperature and many interacting parti- 
cles. 

This section summarizes the two main ideas of this 
paper: the relation between the connectivity of a graph 
and the optical parameter, and the analysis of the power 
of trapping of a completely connected region. 




The Model 
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Let us consider a 2-dimensional rotational invariant 
graph where the connectivity f r<r +i between two layers 
r and r + 1 is not constant (3Cjv case), but decreases as 



fr.r+1 E} l0 pd 



hop^r,r+l ^ 



' N r 



= 1 



N 



(68) 



d rjr +i is the relative degree of a site of the r-th subgraph 
towards the r + Tth subgraph, N r is the number of ver- 
tices in the layer r, and 7 is a parameter that controls 
how fast the connectivity decreases. Note that this choice 
is consistent because we can always tune Eh op to be small 
enough, iV r = TV uniform, and c Lr+ i to be such that we 
can approximate the rhs of Eq. (68). The %n model cor- 
responds to 7 — > 00 together with N. We expect that 
the larger 7 is, the larger is the localization. 



B. Trapped surface 

The rotational symmetry allows us to perform the 
same diagonalization procedure followed in Sec. |IV| That 
is, we will diagonalize the Hamiltonian in blocks corre- 
sponding to the eigenspaces of the rotation operator L. 
Unlike the %n graph, the diagonalization in such sub- 
spaces must be done numerically. This is not a problem 
in the single particle sector. We can then determine the 
ground state and the gap of the system. 

In Fig. [lOj we have found the spectrum of the system. 
We see that the ground state is protected by a gap that 
increases linearly with the size of the completely con- 
nected region N. By the same argument as in the %n 
model, this gap guarantees that the completely connected 
region is a trapped surface robust against finite temper- 
ature and against adding several interacting particles. 

To characterize the ground state of the system, the 
probability distribution of the position of the particle is 
plot in Fig. 11 We note that, unlike the %n model, the 



particle is not completely localized inside the trapped 
surface, but localized in a small region in and around the 
trapping surface. 

In Fig. [I2j the localization of the particle inside the 
completely connected region is plotted against its size 



N 



FIG. 10. Plot of the energies of the ground state Eq, the 
first excited state Ei, and the state with maximum energy 
B m „, with respect to the size of the trapping surface, N, 
for a rotational invariant graph with decaying connectivity 
(1 + N/r), in the single particle sector. Note that the gap 
Ei — Eq increases linearly with N. 




FIG. 11. Probability distribution of the position of the parti- 
cle. The particle is completely confined around the trapping 
surface. The larger the fall-off coefficient 7, the larger the 
confinement. 



and for three different parameters 7 = 1/2, 1 and 2. Both 
this plot and Fig. [IT] illustrate that, the higher the pa- 
rameter 7, the stronger the localization in the completely 
connected region. 



VII. CONCLUSIONS 

Quantum mechanics drastically deviates from the or- 
dinary intuition because of its nonlocality, or action-at- 
distance. Wes do, however, know that in condensed mat- 
ter systems we observe excitations with relativistic dis- 
persion relations. Why is this? An analogue of this prob- 
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FIG. 12. Localization of a single particle in a trapped surface 
vs its size. The localization of the particle increases with the 
size N and the fall-off coefficient 7 (see the definition of 7 in 
Eq. del). 



lem is sound in classical systems. Newtonian mechanics 
is nonlocal but Newton himself evaluated the speed of 
propagation of sound in gases. 

It is now fairly well-understood that theories with 
nonlocal action but local interactions have excitations 
that are locally covariant in certain limits. The amount 
of nonlocality can be bounded using the Lieb-Robinson 
bounds |18) . In this paper we have explored this phe- 
nomenon in the Hubbard model. 

We found that the vertex degree plays a central role as 
it is related to the local speed of propagation of the classi- 
cal probability distribution of matter on the lattice. Our 
analysis is relevant for the idea that relativity, general 
and special, may be emergent from an underlying the- 
ory which is local and quantum. Going further than the 
analysis of the Lieb-Robinson bounds in JTS] , we looked 
at emergence of an effective curved geometry in the Hub- 
bard model. The quantum Hamiltonian evolution of the 
one-particle density for a particular set of states leads to 
an equation of motion that approximates the wave equa- 
tion for a scalar particle in curved space. The same equa- 
tion appears in optical systems with a varying refraction 
index. 

It is widely known that this equation can be cast in 
the domain of general relativity by noticing that it is 
the same as evolution in Gordon's space-time. This con- 
nection cannot be made in 1+1-dimensions, the case we 



considered, but equation (59) can be extended to more 



than one spatial dimension by considering graphs which 
can be foliated in more than one directions. 

Note also that this has been achieved on a Bose- 
Hubbard model on graphs with varying vertex degree and 
multiple edges between sites. The important picture to 
keep in mind to understand this mechanism is Fig.[TJ The 
graph modifies the strength of the interaction of the orig- 
inal Bose-Hubbard model and, for the particular states 



we are interested in, we can obtain Eq. ( 59 1 for the par- 
ticle localization. 

Using similar methods, we confirmed previous argu- 
ments 13J and showed that regions of high connectivity 
in the lattice trap matter. We did this by showing that 
this ground state is protected by a gap which increases 
linearly with the size of the completely connected sub- 
graph. 

As with the emergent wave equation, this result is rem- 
iniscent of the physics of analogue gravity systems . It is 
then interesting to consider the model of (T3] as a rather 
unusal analogue gravity model. 

We will close with the two obvious next questions, 
which we leave for future work: 1) What is the timescale 
of thermalization of matter in the completely connected 
regions? Does it support the conjecture in [T5} [20] that 
black holes are the fastest thermalizers in nature? 2) 
How do our results change if we include the backreaction 
of the matter on the graph? (Our study makes it essen- 
tial in order to have evaporation.) What is the effective 
dynamics of the full Hamiltonian? 
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APPENDIX: EVOLUTION IN THE 
MULTI-PARTICLE SECTOR 

In this Appendix, we generalize the results of Section 
[V]to the case of multiple particles hopping on the graph. 

Free motion. Let us first consider K distinguishable 
particles in states which are tensor products of delocal- 
ized states. This will allow us to track every single par- 
ticle using a number operator. 

Since the particles are distinguishable, we cannot use 
the standard second quantization formalism. Instead, a 
quantum state of the system must be described by 



I* 



N 

mi ,. . . : ttik = 1 



mi <8> 



\m K ) , 



where \rrii) is the state of the i-th particle at the site m^, 
and N is the total number of sites of the lattice. Note 
that the dimension of the Hilbert space is N K . 

Let us also introduce a number operators Ni, defined as 
the tensor product of the number operator on the site of 
the particle i with identity operators on all sites except i 
so that the total number operator is N — 2*=i N%. This 
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number operator counts the total number of particles on 
the graph. To see how to go from the 1-particle case 
to this number operator, let us define Nf as the number 
operator on the vertex k for the particle i. In this case we 
can write N = k TV*, and define a number operator 

for a region A of the graph as Na — J2keAi^i- This 
consistently measures the number of particles in a region 
A of the graph. 

It is easiest to start with the 2-particle case and then 
extend the analysis to K particles. The main complica- 
tion compared to the 1-particle case is that we now have 
a density matrix with 4 indices. In the non-diagonal case, 
we have 

P= E ^fci,fc a ,*i',k a /(*)l*i)<*2|®l*l)(*2l- (69) 

k\ ,k2,ki/,k2' 

Fortunately, we also have H = Hi + H2, with [Hi, H 2 ] = 
0. A straightforward calculation shows that we can use 
the result from the 1-particle case: 



H,[H,p] = Hi,[Hi,p] + H 2l [H 2 ,p] +R, (70) 



where the mixed term is 



R = 2( HiH 2 p + P HiH 2 - HipH 2 - H 2 pH x 



(71) 



We follow the same steps as in the 1-particle case and 
evaluate the second derivative of the expectation value 
of the two number operators to find 



dfal 



dWk 



Tr 



[H u [H llP ]]Nl 



1 

Tr[[H 2 ,[H 2 ,p}}Nl 
Tr \Rp(Nl + N 2 k/ 



(72) 



We can now use the same argument as in the 2-particle 
case to show that all the terms Rij vanish independently. 
We define a k :— ~ ^\ a l k to be the probability of finding 
a particle at k. This quantity satisfies the same equation 
as the one-particle sector probability: 



]Td t V fe (i) = E Vc 4- 



(74) 



That is, particles are independent from each other and 
each follows its own equation ( 74 ) . 



Interaction. Including interaction, in general, the 
Hamiltonian of the system can be written as 



£=E^+E# 



i<3 



where Hi 



(75) 



For what fol- 



ic— i— 1 times 



(74) 



and 

when we add an 



lows, again we first check the case of two particles 
calculate what happens to Eq. 
interaction term. The most genera 
nian is of the form H int 



— 12ijkl Uij 



interaction Hamilto- 
kl a^aj^a^af. We 
make a simplifying assumption which is natural for de- 
localized states, and require that the potential respects 
their symmetries. The simplest such local potential is of 
the form 



H mt =J2h\k)(k\®\k)(k\. 



(76) 



If we add this hamiltonian to the free one, additional 
terms appear on the rhs in Eq. ( 74 ) . These can be traced 



back to the additional commutators in the expansion of 
the full Hamiltonian, 



It is easy to understand what happens when p = p\ ® 
p 2 , for which U kuk2>kl /, k2 ,(t) = U k ^ k? (t)U k ,, k7 ,(t): the 
first and second terms on the r.h.s. of Eq. ( |72| ) reduce 
to discrete second-order derivatives, as in the 1-particle 
case. We therefore only have to deal with the Tr [Rp(N~l + 
N%)] term. 

Another way to see this is by noticing that Nj acts 
as a projector on the one-particle states. We have 
Tr[HipH 2 ] = Tr[H 2 H 1 p], using the properties of the 
trace. Term by term, we can show that, for each Hi 
and Hj , the mixed term vanishes for any number of par- 
ticles. Defining R^ — HiHjp + pHiHj — HipHj — HjpHi, 
we find that the many-particle equation becomes, 



Es t 2 c4W = E v< 4+ 2Tr 



<»i> 



E^ 



(73) 



where Va| is the second order discrete derivative 



d x (cd x (-)). 



C := 



Hi , [H int ,P\ + H int , [Hi , p] 
h(l->2)+ \H m tAHnUp] 



(77) 



that come from expanding the full hamiltonian Ht = 
Hi + H 2 + H^t in the double commutators and keeping 
only the terms involving free Hamiltonians. Since we 
want to distinguish the single particles in the continuum 
limit, we consider the approximation ||-ffmt|| <C H-H1/2II, 



with || TH = sup, 



\\Tv\\ 
\v\ 



This gives a product of the density 



matrices of the single particles. 

A lengthly but straightforward calculation shows that, 
when we restrict to diagonal density matrices in the sub- 
space of the single particles, we find: 

Tr[CNl] = Tv[CN 2 k ] = 0. 

This means that, surprisingly, in the continuum limit, 
the wave equations for the fields are decoupled if we use 



a local potential of the form ( 76 1 
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The next generalization is a potential which is slightly 
non-local. The simplest such potential is given by 

H int = J2h\k + l)(k + l\®\k)(k\ (78) 

k 

for the case of 2-particles. In this case, an even longer but 
still straightforward evaluation of Tr[CN%] or Tr[CW|] 
shows that these traces are nonzero. They take the form 

Ti[CNl] = f k , k U k , k (t)U k+1 , k+1 (t)(I k+1 - J fc ), (79) 
and 

Tr[CN%] = f k , k U k+lik+1 {t)U k!k (t)(I k+1 - J fc ). (80) 

That is, in the continuum limit, the equations for the 
two probability fields are coupled. If we define a k (t) := 
Uk,k{t) and /3fc(t) := U k ,k the probability fields obey 

d?a k - Va fc = a k l3 k+ ifk,k{Ik+i - h) (81) 
df/3 k - V/3 fc = /3 k a k+1 f k , k (I k+1 - h) (82) 

If we define fi(x) to be the continuum equivalent of fk,k 
and I(x) the continuum equivalent of I k , we have 

Ua{x, t) = a(x, t)0(x, t)n{x)I'{x) (83) 
D/3(x, t) = [3{x, t)a(x, t)n{x)I'{x). (84) 

These equations can be straightforwardly generalized to 
the case of more than two particles if the potential is the 



sum of a 2-body interaction for each pair of particles of 
the type H} 3 nt , with 

k »-th 

<g> i <g> • • • <g> i <g> |fc - l) (fc - i| <g>i • • • ® i f k j . 

s v ' 

J-th 

(85) 

The sum (ij) is over all pairs of particles. If a l {x,t) is 
the probability field for particle i, we obtain, for each 
particle, the following coupled differential equations: 

Ua l = a l C^dr 3 a 1 )^. (86) 

We can now define mf = (E j¥i dI ij Q>i)(j,. The EOM 
becomes, 

□a 4 - m?(a;,i)o! < = 

The effective mass mf is determined by the interaction 
with the other particles. The study of the interaction is 
important to understand, but beyond of the scope of the 
paper. 
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